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AN  IMPROVED  CHA  ALGORITHM  FOR  TRACKING  LOW-ANGLE  TARGETS (U) 


B.J.  Rook  and  J.  Litva 


ABSTRACT 

y/^' An  improved  Correlation  Height  Analysis  (CHA)'  algorithm  is 
described  and  further  supporting  evidence  for  its  effectiveness  is 
provided.  The  underlying  principles  are  the  same  as  the  original 
CHA  algorithm,  namely;  (a)  modelling  the  behavior  of  tracking 
radars  in  the  low-angle  region;  (b)  comparing  results  of  the  model 
with  measured  data  in  a  correlation  process ;  and  (c)  selecting  the 
true  track  from  the  family  of  tracks  provided  by  the  correlator. 
The  improvements  to  be  discussed  are  concerned  mainly  with 
increasing  the  efficiency  and  robustness  of  the  algorithm.  The 
most  significant  is  the  procedure  for  effecting  item  (c).  In 
addition,  a  description  is  given  of  the  software  implementation 
of  the  algorithm  together  with  the  sampled-aperture  radar 
facility  used  to  record  data  for  testing  the  algorithm.  CHA 
low-angle  tracking  is  demonstrated  using  experimental  data 
recorded  with  this  tacilityJU) 
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1.  INTRODUCTION 

In  recent  years,  a  considerable  amount  of  research  has  been  devoted  to 
finding  a  solution  to  the  problem  of  tracking  targets  in  the  low-angle  region, 
in  particular,  when  these  targets  are  tracked  over  the  surface  of  the  sea. 

In  this  region,  the  performance  of  conventional  tracking  radars  is  seriously 
degraded  owing  to  the  reflection  of  radar  signals  by  the  surface  of  the  water. 
The  main  thrust  of  the  research  in  this  area  has  been  directed  towards  finding 
methods  for  eliminating  or  at  least  reducing  the  degrading  effect  of  the 
reflected  or  indirect  signal  so  as  to  improve  the  tracking  accuracy  of  the 
radar. 

I.itva1  has  described  a  technique  whereby  advantage  is  taken  of  the 
information  due  to  multipath  contained  in  the  composite  signal  detected  by 
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the  radar.  This  technique,  which  is  known  as  the  Jorrelation  Height  Analysis 
(C1D)  technique,  uses  a  model  to  simulate  the  multipath  propagation  in  this 
region.  The  height  of  the  target  is  derived  by  cross-correlating  theoretical 
data  obtained  from  the  model  with  measured  data.  It  should  be  noted  that  the 
technique  identifies  a  family  of  target  tracks  all  of  which  are  initially 
equally  lkiely. 

In  this  report,  the  formulation  of  an  improved  CHA  algorithm  is 
described.  In  addition,  recent  results  are  presented  to  demonstrate  the 
capability  of  this  algorithm  with  experimental  data. 


2.  ASSUMPTIONS 

In  the  following  mathematical  development  leading  to  the  derivation  of 
the  correlation  function,  simple  ray  theory  is  used  to  describe  the  interfer¬ 
ence  between  two  waves  following  a  direct  and  an  indirect  path.  In  this 
description,  it  is  a=jumed  that  there  is  only  one  reflecting  surface  and  that 
this  surface  is  perpendicular  to  the  vertical  at  the  point  of  reflection.  If, 
in  practice,  it  is  found  that  the  presence  of  surface  tilts  so  degrades  the 
performance  of  the  algorithm,  it  is  further  assumed  that  these  effects  can  be 
made  to  be  negligible  by  assorted  techniques  such  as,  range  averaging  of  the 
radar  signal  and  modelling  the  tilts  using  data  from  sensors  placed  aboard 
the  ship. 


3.  PLANE- EARTH  GEOMETRY 

In  the  development  of  formulae  to  describe  accurately  the  signal 
detected  by  the  radar  in  the  loy-angle  region  we  first  consider  the  plane- 
earth  geometry  which  is  shown  in  Figure  1.  Here,  we  have  a  radar  antenna  at 
height  z±  above  a  horizontal  plane  and  a  target  at  height  Z2  above  the  same 
plane.  The  composite  signal  detected  by  the  radar,  assuming  the  antenna  to 
be  isotropic,  is  the  sum  of  Ej  the  direct  signal  due  to  a  plane-wave  along 
the  direct  path  and  the  indirect  signal  due  to  a  plane-wave  reflected 
from  the  surface  of  the  earth. 


Fiyure  t.  Low-angle  Geometry  (Plane  Farth) 
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For  the  moment,  we  shall  asume  that  the  presence  of  the  plane-earth  ran 
be  represented  by  the  plane-wave  reflection  coefficient  T  which  is  a  function 
of  the  grazing  angle  iju,  and  is,  in  general,  a  complex  quantity  represented- 
by 

r  =  o  e  0  <  P  <  1.  (1) 

The  quantity  o  gives  the  reduction  in  the  amplitude  of  the  signal  on  reflec¬ 
tion  and  o  gives  its  phase  retardation.  In  the  following  sub-section,  the 
composite  signal  observed  at  z^  due  to  the  presence  of  the  direct  and  indirect 
signals  is  formulated  and  discussed. 


3.1  FORMULATION  OF  THE  COMPOSITE  SIGNAL  OBSERVED  AT  z1 


be 


In  formulating  the  composite  signal  at  z^, 
given  by 


E 


d 


j (mt-kr  ) 
E  e  d 


let  the  direct  signal 


(2) 


where 


Eq  =  the  amplitude 

u>  *  the  angular  frequency 
r^  =  the  length  of  the  direct  path,  and 

k  =  is  the  wave  number,  X  being  the  free-space  wavelength. 
Similarly,  let  the  indirect  signal  E^  be  given  as 


.1  Ut-kr  -$) 

E.  =  E  D  e  1  Til 

l  o 


where  r^  is  the  length  of  the  indirect  path.  The  composite  signal  observed 
at  z^  is  simply  the  sum  of  E^  and  E^,  thus 
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o  e 


-j  { <f-*-k(ri-rd)  > 


(4) 


Equation  (4)  defines  an  interference  pattern  in  space  whose  shape  is 
given  by  the  quantity  inside  the  square  brackets.  The  factor  outside  the 
brackets  gives  the  signal  that  would  be  observed  if  only  the  direct  ray  were 
present.  This  interference  pattern  is  of  considerable  importance  in  that  its 
detailed  shape  depends  upon  the  location  of  the  target.  Insofar  as  modelling 
the  multipath  propagation  in  the  low-angle  region  (interference  region)  is 
concerned,  our  only  interest  is  in  the  amplitude  and  phase  of  the  indirect 
signal  relative  to  the  direct  signal,  consequently,  the  factor  outside  the 

£  ^  , 
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square  brackets  may  be  discarded.  Therefore,  eqn.  (4)  ma*  be  rewritten 
g-ving 


-j{$+k(r  -r  ) ) 

E  -  1  +  o  e  1  “  (5) 

Thus,  eqn.  (5)  describes,  in  complex  fcrm,  the  interference  pattern  at  z^  due 
to  the  superposition  of  two  signals  which  have  related  amplitude  ai  i  phase 
characteristics. 

From  the  geometry  of  Figure  1,  tie  lengths  of  the  direct  and  indirect 
paths  may  be  determined  to  be 

rd  "  (r2  +  <V2i)2]1/2  (6) 

and 

r£  -  [r2  +  (z1+z2)2]1/r2  ,  (7; 

therefore,  the  path-length  difference  Ar  »  r^-r^  may  be  determined  by  expand¬ 
ing  each  path  length  by  the  binomial  theorem  and  subtracting,  to  wit 


Ar  * 


2z1Z2 


1 


z2xz2 

zrz2 

2r2 


(8) 


This  expression  is  unwieldly,  but  fortunately  in  the  low-angle  region,  that 
is,  where  <l>i  and  are  snail,  the  expression  may  be  simplified  to  the  more 
usual  form  of 


Ar  as 


2z  z 
12 


(9) 


with  r  =  r^,  the  slant  or  measured  radar  range  to  the  target.  Therefore,  we 
may  rewrite  eqn.  (5) ,  giving 


E  -  1  +  p  e-J(*+Wr)  (10) 

Thus,  eqns.  (9)  and  (10)  map  out  the  phase  and  amplitude  of  an  inter¬ 
ference  pattern  observed  ac  z^  as  a  function  of  r  and  Z2-  Alternatively,  if 
r  and  Z2  are  considered  to  be  constant,  eqns.  (9)  and  (10)  map  out  the  phase 
and  amplitude  of  an  interference  pattern  as  a  function  of  z^.  From  this  it 
follows  that  if  a  number  of  elements  were  vertically  disposed  above  the 
horizontal  plane  at  the  receiving  end,  we  may,  at  least  conceptually,  sample 
simultaneously  the  phase  and  amplitude  at  the  element  positions  for  given 
values  of  r  and  Z2»  Thus,  we  have  introduced  the  concept  of  the  sampled 
aperti  re,  whereby  we  sample  the  composite  signal  at  discrete  intervals  on  a 
finite  straight  line  in  space  at  some  instant  of  time  where,  at  this  instant, 
r  and  Z2  are  considered  to  be  constant.  The  concept  of  the  sampled  aperture 
and  its  application  to  solving  the  low-angle  tracking  problem  is  of  consider¬ 
able  importance,  and  will  be  discussed  in  much  greater  detail  later  in  this 
report. 
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Suffice  it  to  say  that  in  order  to  describe  the  relative  phase  and 
amplitude  distribution  across  this  aperture,  eqn.  (10)  must  be  evaluated. 
However,  this  equation  cannot  be  evaluated  until  a  complete  knowledge  of  the 
behavior  ot  p  and  is  known.  In  addition,  since  a  sampled  aperture  is  being 
considered,  the  distribution  of  phase  across  this  aperture  for  varying  angles- 
of-arrival  of  the  direct  and  indirect  signals  must  also  be  known. 


3.2  DETERMINATION  OF  THE  ANGLES-OF-ARRIVAL  OF  THE  DIRECT  AND  INDIRECT  SIGNALS 


To  fully  describe  the  resultant  phase  and  amplitude  distribution  across 
the  sampled  aperture  for  a  given  r  and  Z2,  we  must  be  able  to  calculate  the 
angles-of-arrival  of  the  direct  and  indirect  signals  relative  to  the  boresight 
of  the  aperture.  After  carrying  out  these  calculations,  we  then  find  that  we 
are  able  to  derive  the  quantities  p  and  $  which,  in  turn,  leads  to  an  evalua¬ 
tion  of  eqn.  (10). 


In  deriving  formulae  to  determine  the  angles-of-arrival  of  the  direct 
signal  ”j  and  the  indirect  signal  9^  relative  to  the  boresight  of  the  apertur-; 
it  is  necessary  to  establish  a  convention  for  signs  for  the  angles  shown  in 
Figure  1.  The  angles  pj,  *2  and  9p  (the  pointing  angle)  are  measured  positive 
upwards  from  the  horizontal  line,  whereas  the  angles  84  and  9j  are  measured 
positive  upwards  from  the  aperture  boresight.  In  the  construction  of  Figure 
1,  both  6^  and  9^  arc  shown  to  be  numerically  negative.  It  follows  that 


and 


tanik^ 


tan«<2 


VZ2 


(ID 


(12) 


In  the  low-angle  region,  where  and  ^  are  small,  eqn.  (11)  and  eqn.  (12) 
reduce  to 


and 


Z2_Z1 


v2  r 


(13) 


(14) 


whereupon  the  angles-of-arrival  of  the  direct  si^-al 
signal  relative  to  the  aperture  boresight,  may  be 


9,  and  the  indirect 
given  as 


and 


ik. -9 
1  P 


(15) 


■(<k_+9  ) 
2  P 


(16) 


From  the  preceding  we  note  that  in  the  process  of  deriving  9.  and  9.  we  have 
also  obtained  the  grazing  angle 
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3.3  DETERMINATION  OF  THE  PLANE-EARTH  REFLECTION  COEFFICIENT 

The  plane-earth  reflection  coefficient,  which  is  in  general  a  complex 
quantity,  is  given  by  eqn.  (1).  The  evaluation  of  this  equation  is  based 
primarily  on  a  knowledge  of  the  grazing  angle  <i>2  and  the  electrical  properties 
of  the  reflecting  surface.  For  a  smooth  plane  surface,  the  complex  reflection 
coefficient  may  be  given  for  either  vertically  or  horizontally  polarized  waves 
as 


o  e 
v 


-J*v  ccsin(n»2)  -  -cos^(u>2) 

£csir.(ii2)  +^ec- cos2(y2) 


(17) 


and 


•H  H  PH  6 


-J$H  ain(*2)  - cos2(*2) 
sin(W2>  +  ^ec  -  cos^*  ) 


(18) 


where  the  subscripts  V  and  H  denote  vertical  and  horizontal  polarization 
respectively. 

The  term  ec  in  eqns.  (17)— (18)  is  known  as  the  complex  dielectric  constant  of 
the  surface  and  is  given  as 


e  -  -  j  -2_ 

c  €o 


(19) 


where  e, cq  »  the  ratio  of  the  dielectric  constant  of  the  surface  to  that  of 
free  space 

o  »  the  conductivity  of  the  surface  in  mho/m 
- 12 

tQ  “  8.84  x  10  Farads/m  is  the  dielectric  constant  of  free  space. 
Alternatively,  eqn.  (19)  may  be  expressed  in  terms  of  the  wavelength  X  as 


e  * - j  60Xo 

c  e 

o 


(19a) 


provided  \  is  expressed  in  metres. 

The  values  of  l/.;q  and  a  to  be  used  in  eqn.  (19)  and  (19a)  are  determined  by 
the  characteristics  of  the  surface  that  is  under  consideration.  Values  for 
r/e0  and  c  over  sea-water  and  fresh-water  surfaces  are  given  in  Table  1.  The 
electrical  properties  of  these  two  surfaces  were  chosen  since  the-  . j:  -  .sent 
the  extremes  of  the  range  of  water  surfaces  that  are  likely  to  be  i  ,ci  .  ?red 
in  practice. 
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TABLE  1 

Electrical  Properties  of  Water 


Surface 


Sea-water 


Fresh-water 


T~'  11 - 

0 

mho/m 

cle0 

4  3 

80 

10‘2 

80 

If  we  note  that  in  the  low-angle  region  ^2  *s  small  and  that  for  either 
sea-water  or  fresh-water  jecJ  >>  i,  then  eqn.  (17)  and  (18)  may  be  simplified 
to 


'j0v  a  Vi'fc-  1 
+  1 


H  *  °H  6 


'*'2  ~\Cc 

*2  +yTc 


Figure  2  shows  plots  of  the  magnitude  of  the  reflection  coefficient  0 
against  grazing  angle  V2  for  vertically  and  horizontally  polarized  waves  over 
a  perfectly  smooth  sea-surface.  Figure  3  gives  plots  of  the  phase  of  the 
reflection  coefficient  for  vertically  and  horizontally  polarized  waves.  These 
results  were  calculated  with  eqns.  (17)— (18)  using  a  wavelength  X  *  0.032m 

(X— hand) . 


4.  CURVED-EARTH  GEOMETRY 

As  the  value  of  r  increases,  the  foregoing  mathematical  expressions 
become  increasingly  inaccurate  because  of  the  curvature  of  the  earth.  In 
order  to  correct  for  this  curvature,  we  must  now  consider  the  curved-earth 
geometry  given  in  Figure  4.  Insofar  as  the  low-angle  region  (interference 
region)  is  concerned,  the  diagram  of  Figure  4  is  grossly  exagerated, 
however,  it  has  been  presented  in  this  manner  in  order  to  clearly  define 
the  quantities  of  interest. 

Since  we  are  dealing  with  a  curved-earth  rather  than  a  plane-earth 
geometry  a  new  quantity,  1  amely,  the  divergence  factor  must  be  introduced. 

The  divergence  factor,  denoted  D,  gives  a  measure  of  the  defocussing  that 
occurs  when  an  electromegnetic  wave  is  reflected  from  a  curved  surface.  How- 
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smooth  sea 
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Figt.  3.  Verticil  and  Horizontal  Reflection  Coefficient  Phases 


1^  - 


ever,  before  this  an-4  ocher  related  quantities  may  be  determined,  the  distance 
from  either  terminal  point  to  the  reflecting  point  must  first  be  calculated. 


4.1  DERIVATION  OF  THE  POINT  OF  REFLECTION  ON  A  CURVED  SURFACE 

When  dealing  with  a  plane-earth  geometry  we  found  that  the  parameters 
Chat  had  the  greatest  impact  when  deriving  interference  patterns  (see  eqn. 
(10))  were  the  path-length  difference  Ar  and  the  grazing  angle  This  is 

also  true  for  a  curved-earth  geometry,  but  first  we  must  determine  the 
position  of  the  reflecting  point,  or  how  r  is  broken  into  the  length  rj  and 
r 2-  Once  r^  and  T2  have  been  determined,  corrected  values  for  Ar  and  ^  may 
be  calculated. 

If  we  are  given  the  values  of  z ^  and  r,  it  is  possible  to  find  r^ 
by  solving  the  cubic  equation4 

f (r)  -  2r?-3rjr  +  lr;  -  2A  (z.+z_)]r.  +  2A  z.r  -  0,  (22) 

11  eid£  i  ei 

where  A  is  the  effective  earth's  radius,  normally  taken  to  be  8.5  x  I0&m. 
e 

This  cubic  equation  has  the  formal  solution: 

ri  - 1 +  p  cos  {—)  • 


where 


10 


f 


r 


S3  1  e 


(*1+Z2)  + 


(!) 


(23) 


and 


«  » 


Unfortunately,  the  formal  solution  to  this  equation  is  not  very 
practical  for  implementation  on  a  digital  computer,  particularly  when  one 
considers  that  there  is  an  inverse  cosine  to  be  evaluated.  A  more  practical 
method  for  its  solution  is  described  in  the  following  suh- section. 


4.2  NEWTON1 S  METHOD  FOR  SOLVING  THE  CUBIC  EQUATION 

An  efficient  method  for  solving  the  cubic  equation  with  a  digital 
computer  is  by  using  Newton's5  method.  When  written  in  its  recursive  form, 
Newton's  method  take:  the  form 


f[(OJ 

^1^+1  “  (rl5n  ~  f '  [  (r  )  ]  * 

1  n 


(24) 


The  parentheses  surrounding  r^  are  to  distinguish  the  subscript  of  r  from 
the  index  n  of  the  recursion  formula. 


Initially,  (ri)n„^  is  assigned  a  value  which  is  near  the  desired  root 
of  the  cubic  equation.  This  value  is  determined  by  noting  from  the  olane- 
earth  geometry  of  Figure  1  that 


^n-l 


(25) 


Equation  (24)  is  then  computed  in  an  iterative  manner  until  | (r. )  -(r^)n+1 
<_  1.0m.  In  practice,  three  iterations  of  eqn.  (24)  are  required  for 
this  inequality  to  be  satisfied. 


4.3  DETERMINATION  OF  THE  PATH-LENGTH  DIFFERENCE 

Once  the  position  of  the  reflecting  point  has  been  determined,  we  may 
then  compute  all  other  quantities  of  interest  relative  to  the  plane  which  is 
tangent  to  the  surface  of  the  earth  ar  r^. 

To  begin,  we  calculate  new  values  of  z^  and  zj*  namely  z[  and  z^. 
Referring  to  Figure  4,  it  is  convenient  to  define  parameters  li  and  8.2  which 
give  distances  from  the  tangent  plane  to  the  earth's  surface.  For  small 
angles  they  are 


r  s* 
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and 


(26) 

(27) 


In  the  low-angle  region  where  and  Vo  are  small,  the  distance  along  the 
tangent  plane  from  the  respective  terminal  points  to  the  reflection  point  may 
be  taken  to  be  r^  and  r2-  In  a  straight  forward  manner*'  it  can  be  shown  that 


1,2  2A 


from  which  the  new  heights  z^  and  z^  may  be  determined  thus, 


z  —  z  —  — 
1  1  2A 


and 


<  2 
Z  f  —  ■ 

2  "  2  2A 


(28) 


1.29) 


(30 


Be  referring  to  Figures  (1)  and  (4)  and  using  eqn.  (9),  the  corrected  path- 
length  difference  can  be  written  as 


hr  2 


-  »  • 

2Z1Z2 


(31) 


4.4  DETERMINATION  OF  6.  and  9. 

d  l 

The  calculation  of  the  angles-of-ar  '  val  for  the  direct  and  indirect 
signals,  6^  and  9^  respectively,  are  carried  out  in  a  manner  similar  to  that 
for  the  plane-earth  case,  with  the  exception  of  now  having  to  consider  n, 
which  is  the  angle  between  the  horizontal  plane  and  a  plane  parallel  to  the 
tangent  plane. 

It  follows  from  an  inspection  of  Figure  4  that 


n 


Furthermore, 


1 


Z2_Z1 


r 


(32) 


(33) 


and 


12 


(34) 


whereupon  the  angles-of-arrival  of  the  direct  and  the  indirect  signals  relative 
to  the  aperture  boresight  are 


and 


0 


d 


4,^-6 


P 


n 


o  -  +  0  +  n). 

i  i  p 


(35) 


(36) 


4.5  THE  DIVERGENCE  FACTOR 


When  a  bundle  of  parallel  rays  is  incident  on  the  convex  side  of  a 
curved  surface,  the  rays  diverge  after  reflection  because  each  has  a  slightly 
different  angle  of  incidence.  Because  of  this  divergence,  the  reflected 
signal  is  defocussed  and  the  power  density  is  reduced.  The  divergence  factor 
can  be  derived  solely  from  geometrical  considerations7  and  is  given  here  in  a 
form  which  incorporates  the  usual  approximations  which  are  applicable  in  the 
low-angle  region,  namely 


D  s 


1 


2rif2 

V*2 


(37) 


5.  SURFACE  ROUGHNESS  FACTOR 

The  discussion  thus  far  has  been  based  on  the  assumption  that  the 
surface  was  smooth  in  the  neighborhood  of  the  specular  reflecting  point.  We 
now  consider  the  modifications  brought  about  by  the  introduction  of  a  rough 
surface. 

In  general,  when  parallel  rays  are  reflected  by  a  rough  surface,  they 
undergo  changes  in  path  length  and  amplitude  because  of  the  deviations  from 
a  smooth  reflecting  surface.  These  changes  result  in  a  portion  of  the 
reflected  signal  becoming  diffuse.  Where  previously  it  was  totally  coherent 
it  now  consists  of  two  components:  (1)  a  diffuse  component  and  (2)  a  coherent 
component  with  reduced  magnitude. 

The  reduction  to  the  magnitude  of  the  coherent  component  brought  about 
by  reflections  from  a  rough  surface  is  related  to  the  angle  of  incidence  or 
grazing  angle  and  the  radar  wavelength.  It  is  best  expressed  in  terms  of 
the  RMS  specular  scattering  coefficient  S  which  is  given  by8 

(  j2rausint|;J2) 


S  ■  exp 


(38) 
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where 


the  rms  surface-height 


variation 


(assumed  normally  distr ibuted) . 


For  smaLl  values  of  i^2.  it  is  customary  to  replace  sinv'2  by  -2  and  P’ot 
3  as  a  function  of  the  'apparent  roughness  factor*  A  plot  of  S 

against  0^2/  '  expressed  in  milliradians  is  shown  in  Figure  5. 


Board '  has  recorded  data  that  shows  good  agreement  with  result-, 
predicted  by  eqn.  (38)  for  values  of  0^2/*  1  100  milliradians.  However, 
beyond  100  milliradians,  his  results,  which  are  shown  in  Figure  6*,  imicaies 
that  S  decreases  more  slowly  than  predicted  by  theory. 


Thus,  based  on  the  results  of  Beard,  the  specular  scattering  coefficient 
given  by  eqn.  (38)  is  used  to  determine  the  magnitude  of  the  coherent 
component  of  the  reflected  signal  for  values  of  ^  100  milliradians. 

Beyond  this  value  the  approximate  expression 

3  2  expj-[0. 165(0^2^)  ‘  +  7*^2  +  0.0468]|  (39) 


is  used.  A  plot  of  S  against  aH*s/x*  using  eqns.  (38)  and  (39),  is  shown  in 
Figure  7. 


*  Beard  uses  the  symbols  C/Dr.  h  and  1 h  for  the  parameters  identified  by  S.  oH  and  ^2  respectively  in  this 
report. 
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Figure  5.  Specular  Scattering  Coefficient 
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5.1  RELATIONSHIP  BETWEEN  SEA-STATE  AND 

The  parameter  used  for  quantifying  the  roughness  of  an  ocean  surface  is 
called  'sea-state'.  It  is  related  to  the  rms  surface  height  variation  -ty 
through  the  'significant  wave  height*  H.  In  Nathanson10  the  significant  wave 
height  is  defined  as  the  average  of  the  highest  third  of  the  waves  in  a  wave 
group  and  is  shown  to  be  related  to  ay  by  the  following: 

H  =  4oH  .  (40; 

Thus,  if  we  have  an  estimate  of  the  sea-state  we  can,  by  using  eqn.  (40)  and 
Figure  7,  derive  the  magnitude  of  the  coherent  component  of  the  reflected 
signal  for  varying  radar-target  geometries.  In  Table  2  is  given  a  list  of 
sea-state  numbers  (SS)  with  corresponding  values  of  significant  wave 
height  (H). 


TABLE  2 


Relationship  Between  Sea-state  and  Significant  Wave  Height 


Description  of  Sea 

Sea-State  (SS)  Number 

H  Metres 

Calm 

0 

0.0 

Smooth 

1 

0.1 

Slight 

2 

0.4 

Moderate 

3 

0.8 

Rough 

4 

1.6 

Very  Rough 

5 

3.0 

High 

6 

4.9 

6.  DERIVATION  OF  THE  SIGNAL  AT  Z]  FOR  A  CURVED  ROUGH  SURFACE 

In  the  discussion  in  Sections  4  and  5  a  number  of  new  parameters  were 
introduced  to  describe  a  curved-earth  geometry  and  reflection  from  a  rough 
surface.  These  quantities  are  as  follows: 

(1)  Yr  s  2z^z^/r  -  curved-earth  path-length  difference. 

(2)  D  -  divergence  factor,  to  describe  defocussing  by  a  curved 
surface. 

(3)  S  -  specular  scattering  coefficient,  to  quantify  the  coherency  of 
the  signal  reflected  by  a  rough  surface. 

(4)  Soe  -  reflection  coefficient,  as  modified  by  a  rough  surface. 

By  direct  substitution  into  eqn.  (10)  using  the  items  listed  above  it 
follows  that: 


% 


t  ■ 


RELATIVE  AMPLITUDE  AND  PHASE  (RADIANS) 
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where 


1  +  pDS  e 


a  »  $+kAr. 


Equation  (41)  describes  the  interference  pattern  at  z^  resulting  from  the 
superposition  of, the  direct  signal  and  the  coherent  component  of  the 
reflected  signal.  In  its  description  it  includes  the  effects  of  a  curved- 
earth  geometry  and  reflections  from  a  rough  sea-surface. 

Figure  8  gives  a  plot  of  the  amplitude  and  phase  of  the  interference 
pattern  (E)  as  a  function  of  for  Z2“4.8m,  SS“0  and  r“7000m.  It  should  be 
noted  that  the  interference  pattern  is  periodic  in  z 2  with  a  period  of  about 
26.7m.  Another  interference  pattern  is  given  in  Figure  9  -  this  time  plotted 
versus  range  for  the  fixed  conditions  z^»4.8m,  Z2“10m  and  SS«0.  Both  Figures 
8  and  9  were  derived  using  eqn.  (41).  Rough-sea  results  would  be  similar  to 
these  except  for  a  reduced  amplitude  and  a  reduction  in  the  phase  deviation 
due  to  the  effect  of  the  factor  S  in  eqn.  (41). 


7.  THE  SAMPLED  APERTURE 

The  sampled  aperture  concept,  as  implied  by  its  name,  entails  sampling 
the  amplitude  and  phase  of  a  complex  wavefront  at  N  discrete  points  along  a 
finite  straight  line  in  space  (antenna  aperture).  Using  two-dimensional 
Fourier  processing  it  is  possible  to  develop  N  Independent  beams  and  track 
up  to  N  targets.  The  far  ranging  capabilities  of  sampled- aperture  techniques 
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Figure  8.  Interference  Pattern  Phase  and  Amplitude  for  a  Smooth  See;  Zj  =  4.8m,  r  =  7000.0m. 
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Figure  9.  Interference  Pattern  Phase  and  Amplitude  for  a  Smooth  Sea;  =  4.8m.  =  10.0m. 

form  a  compelling  reason  for  applying  this  technology  to  low-angle  tracking. 
Although  a  full  description  of  possible  applications  is  beyond  the  scope  of 
this  report,  expressions  will  be  developed  for  limited  applications.  An 
example  of  the  latter  is  the  use  of  a  sampled-aperture  radar  to  simulate  a 
standard  fire-control  radar,  a  use  which  might  be  considered  to  occupy  the 
lowest  rung  of  possible  applications. 

Consider  the  linear  array  shown  in  Figure  10  consisting  of  N  equally- 
spaced  elements  separated  by  d.  It  is  convenient  to  consider  arrays  consis¬ 
ting  of  odd  numbers  of  elements;  with  an  odd  number,  the  centre  of  the  array 
coincides  with  an  antenna  element  which  is  labelled  n=*0.  With  this  arrange¬ 
ment,  there  are  (N-l)/2  equal  elements  on  either  side  of  the  centre  of 
reference  element.  If  a  signal  arrives  at  angle  8  relative  to  the  array 
boresight,  then  the  electric  field  distribution  along  the  array  is 

E(n)  -  eJ (ut+kndsin8)  (43) 

where  the  amplitude  is  assumed  to  be  unity. 

Consider  now  that  the  boresight  of  the  array  is  parallel  to  the 
horizontal  plane  in  Figure  4  (i.e.,  9  :0) ,  recalling  that  the  uppermost 
element  is  n  -  (N-l)/2,  ana  that  the  direct  and  indirect  signals  arrive  at 
angles  9^  and  9^  respectively,  then  the  amplitude  and  phase  distribution  of 
the  composite  signal  across  the  array  is 

[  jkndsin9  jkndsin9  1 

E(n)  -  E^(n)+E^(n)  -  le  +  e  je3  ,  (44) 
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Figure  10.  Sampled  Aperture  Geometry 


where  both  signal  amplitudes  are  set  equal  to  one  and  the  phase-path 
difference  between  the  direct  and  indirect  signals  at  the  reference  element 
is  assumed  equal  to  zero.  As  before,  only  that  portion  of  eqn.  (44)  within 
the  square  brackets  is  of  interest  to  us,  therefore,  we  retain  this  portion 
and  discard  the  remainder,  giving 


jkndsinO 


jkndsin8 


We  now  wish  to  incorporate  into  eqn.  (45)  the  parameters,  o,  D,  S  and  a 
which,  for  a  given  r  and  Z2>  are  derived  soley  for  a  value  of  z,  which 
corresponds  to  the  reference  element  (n»0).  For  all  other  n,  the  computed 
values  of  p,  D,  S  and  a  may  be  regarded  as  being  approximately  equal  to  those 
computed  for  n»0  since  the  separation  of  the  elements  of  the  array  is  small 
compared  to  r^  in  Figure  (4).  When  eqn.  (45)  is  rewritten  to  incorporate  the 
above  parameters  and  to  reflect  its  dependence  on  n,  r,  and  Z2,  it  takes  the 


E(n,r,z2)  2  e 


where  for  convenience 


jkndsine . 


+  A(o,r,z2)e 


-j  [0(0, r ,z2)-kndsin9iJ 


A(o,r,z2)  -  p(o,r,z2)*D(o,r,z2)»S(o,r,z2), 


It 
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Alternatively,  eqn.  (46)  may  be  expressed  in  polar  form  with  ampli'ude  and 
phase  of  the  resultant  given  respectively  by 

|E(n,r,22)|  s  jl  +  2A(o,r,22)costknd(sin6d-sin61>  +  n(o,r,z2)] 


1 1/2 


+  A2(o,r,z_) 


and 


$(n,r,z2)  a  tan 


-1 


sin(kndsinOd)-A(o,r ,z2)sin[a(o,r ,z2)-kndsin6i] 


cos(kndsin0  .)+A(o,r,z„)cos[a(o,r ,z_)-kndsin6 , ] 
a  l  l  I 


(Ah) 


(49) 


When  n-0, 


E(o,r,z2)|  -  jl  +  2A(o,r,z2)cos[a(o,r,z2)  ]  +  A2(o,rtz2>j 


1/2 


(50) 


and 


4>(o,r,22)  -  tan 


P  -A(o,r,z2)s:  n[a(o,r,z2)]  "j 
I  1  +  A(o,r,z2)cos[a(o,r,z2)] I* 


(51) 


It  follows  from  eqns.  (48)-(51)  that  the  amplitude  and  phase  distribution 
across  the  aperture  after  normalization  by  the  amplitude  and  phase  at  the 
reference  element  is 


F(n,r,z2) 


E(n,r,z2) 

E(o,r,z2> 


^»(n,r>z2)-i|i(o,r>z2). 


(52) 


It  will  be  noted  in  eqn.  (52)  that  when  n*0,  F(n,r,Z2)  reduces  to  unity  and 
the  corresponding  phase  is  zero.  This  equation  is  important  because  it  gives 
the  resultant  normalized  phase  and  amplitude  distribution  of  a  composite 
wavefront  whose  components  consist  of  the  direct  and  indirect  signals 
arriving  from  a  low-angle  target  located  at  (r,Z2). 

In  the  development  of  formulae  to  describe  the  relative  phase  and 
amplitude  observed  at  z^  due  to  two  signals  (see  eqn.  (4)),  it  was  only 
necessary  to  consider  the  quantity  which  described  the  interference  pattern. 
Thus,  the  remaining  factor  was  discarded  for  ease  of  manipulation  and  clarity 
in  subsequent  mathematical  development.  In  any  event,  if  this  factor  had 
been  included  throughout  the  development,  it  would  have  finally  disappeared 
at  this  noint  due  to  the  normalization  procedure  described  by  eqn.  (52).  In 
additic  this  procedure  effectively  reduces  the  dependency  of  F(n,r,Z2)  on 
range  target  cross-section. 


.  choice  of  the  reference  element  in  the  aperture  is  completely 
ar...  i  .ry  subject  to  the  constraint  of  there  being  an  equal  number  of 

f  4 
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,-lements  on  either  side  of  the  reference.  One  can,  for  example,  divide  the 
array  into  smaller  segments  consisting  of  odd  numbers  of  elements.  This 
allows  sampling  the  wavefront  for  different  values  of  z^,  where  z^  is  the 
height  of  the  centre  element  of  the  array  segment.  It  will  be  shown  in 
Section  8  that  this  flexibility  can  be  used  as  a  basis  for  a  technique  for 
resolving  target  height  ambiguities  which  result  when  low-angle  data  are 
processed  by  the  CHA  filter. 

Although  the  sampled-aperture  technique  has  been  discussed  for  a 
specific  application  i.e.,  the  low-angle  tracking  problem,  one  is  at 
liberty  to  process  the  data  in  a  variety  of  ways.  For  example,  beamforming 
studies  could  be  carried  out  using  the  discrete  Fourier  transform11.  White’s1 
double  null  technique  or  adaptive  array  processing13.  It  is  thought  that  the 
possibilities  for  processing  sampled-aperture  data  are  virtually  boundless 
and  probably  only  limited  by  one's  imagination. 


8.  COMPLEX  CORRELATION  COEFFICIENT 

When  tracking  a  low-angle  target  with  a  sampled-aperture  radar,  the 
target  range  and  antenna  height  can  be  considered  to  be  known  parameters. 
The  unknown  is  z~  and  it  is  our  intent  to  solve  for  z 2  by  cross-correlating 
experimental  data  with  simulated  data. 

Let  us  define  the  measured  value  of  the  signal  at  element  n  to  be 
W(n,r).  Its  normalized  value  T(n,r)  is  then  given  by 


T(n,r) 


w(o,r) 


L± 


(53) 


where  |W(n,r)|  and  $(n,r)  are  the  values  of  amplitude  and  phase  measured  at 
element  n.  When  expressed  in  complex  notation,  eqn.  (53)  becomes 


T(n,r) 


W(n,r) j  ji*(n,r)-0(o,r) 
W(o , r) I  *  l 


(5^) 


For  a  given  value  of  r,  simulated  data  may  be  calculated  for  various  values 
of  Zj  from  eqn.  (52).  Expressed  in  complex  notation,  eqn.  (52)  becomes 


j($(n,r,z  )-0(o,r,z  )} 

«  .  (55) 


The  technique  used  in  solving  for  Z2  consists  of  deriving  cross- 
correlation  coefficients  for  the  measured  amplitude  and  phase  distribution 
across  the  array  aperture,  for  a  particular  r,  with  simulated  data  calculated 
for  closely  spaced  values  of  z~.  Peaks  in  the  cross-correlation  coefficient 
are  used  to  define  possible  values  for  the  target  height  Z2.  Typically,  Z2 
is  sampled  at  1m  intervals  for  a  total  of  about  100  samples.  After  the 
measured  and  simulated  data  are  correlated  the  number  of  possible  target 
heights  may  be  reduced  to  10. 


F(n,r,z2)  - 


E(n,r,z2) 


E(o,r,z2) 


I 


« 
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Using  the  general  expression  given  in  Litva1  the  complex  cor 
coefficient,  denoted  by  C,  roay  be  written  in  terms  of  r  and  as 

(N-l)/2 

C(r,z  )  =■  :  T(n,r)F*(n,r,z  )/o(r)*o(r,z  ), 

1  -(N-l)/2 


(5e>) 


where 


o(r) 


and 


o<r,z7)  * 

and  *  denotes  the  complex  conjugate.  It  should  be  noted  that  the  complex 
correlation  coefficient  given  by  eqn.  (56)  is  calculated  over  all  N  spatially 
sampled  points,  where  N  is  the  total  number  of  antenna  elements  or  a  subset 
of  that  number  in  the  case  where  the  array  is  segmented  into  sub-arrays. 

In  the  following  sub-section,  examples  of  the  complex  correlation 
coefficient,  using  simulated  data,  are  presented  and  examined. 
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T(n,r)T*(n,r) 


(57) 
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F(n,r,z2)F*(n,r,z2)| 


(5c 


8.1  DISCUSSION  OF  SIMULATED  RESULTS 

The  results  to  be  presented  in  this  section  were  derived  using 
simulated  data  for  a  sampled  aperture  consisting  of  seven  elements  with  an 
inter-element  spacing  d  of  0.128m  (Al).  The  boresight  of  the  aperture  is 
assumed  to  be  parallel  to  the  horizontal  plane  (i.e.,  t)  ^0)  with  Z2=10m, 
r=7000m  and  SS»0.  At  thi3  stage  of  the  development  of  the  CHA  algorithm 
only  the  real  part  of  the  complex  correlation  coefficient  is  considered. 
Further  studies  involving  the  imaginary  part  will  be  carried  out  at  a  later 
date. 


Figure  11  gives  the  real  part  of  the  comp'  -.x  correction  coefficient, 
R[C(r,zD],  as  a  function  of  This  result  was  derived  using  all  seven 

elements  of  the  sampled  aperture,  whereby  N»7  and  z^-»4.8m,  where,  it  will 
be  recalled,  z^  gives  the  height  of  the  reference  element. 

It  is  noted  that  R[C(r,Z2)J  repeats  periodically  with  increasing  Z2, 
which  is  consistent  with  the  behavior  of  the  interference  pattern  given  in 
Figure  8.  In  addition,  the  maxima  of  the  correlation  peaks  gradually 
decrease  with  increasing  values  of  tne  index  K,  where  K  identifies  the 
correlatio  i  peaks.  The  true  peak  is  given  by  K»0  and  the  false  or  alias 
peaks  by  K>0.  In  Figure  ll,  one  can  see  the  difficulty  in  distinguishing 
the  true  peak  from  the  adjacent  peaks  simply  on  the  basis  of  its  size. 


Z2  (METRES) 

Figure  1 1.  Synthetic  Cross-correlation  coefficient ;  Zj  =  4.8m.  Z2  =  10.0m.  N  =  7.  0o  =  O.Omr.  r  =  7000.0m. 

In  general,  the  size  of  the  correlation  peaks  does  not  form  a  particu¬ 
larly  goo-  criterion  f >r  selecting  the  true  peak.  The  rate  of  decrease  of 
the  maxima  for  increas  ,ng  K  relies  heavily  or.  the  ratio  of  the  aperture 
length  to  the  scale  length  of  the  interference  pattern.  This  point  is 
brought  out  in  Figure  12,  where  only  the  uppermost  three  elements  (N=3)  are 
used  i:i  deriving  the  complex  correlation  coefficient.  For  this  result  zj  = 
b.056m.  One  sees  that  the  fall-off  in  the  maxima  of  the  peaks  is  much  more 
gradual  and  it  has  become  more  difficult  to  distinguish  the  peak  having  the 
largest  value. 

Whereas,  selection  of  the  true  peak  on  the  basis  of  its  maximum  value 
is  unreliable,  the  magnitudes  of  the  peaks  can  be  used  as  a  criterion  to 
reduce  the  number  of  alias  tracks.  A  reference  level  can  be  set  and  peaks 
can  be  excluded  on  the  basis  of  their  maxima  falling  below  the  reference 
level.  As  an  example  of  this,  Figure  13  shows  four  target  tracks,  the  true 
track  identified  by  K»0  and  three  alias  tracks  identified  by  K**l ,  K*2,  and 
K=3;  these  tracks  were  selected  on  the  basis  of  the  maxima  of  their  correla¬ 
tion  peaks  exceeding  0.8.  All  remaining  alias  trackn  have  maxima  that  fall 
short  of  0.8  and  are  rejected  on  that  basis. 

It  is  noted  that  the  alias  tracks  converge  to  the  true  track,  i.e. , 
Z2=*l0m,  at  r”0m.  This  behavior  illustrates  the  decrease  in  the  scale  length 
of  the  interference  pattern  for  decreasing  range.  It  follows  that  if  r  is 
decreased  in  Figure  11,  then  the  peaks  will  crowd  together. 

These  results  corroborate  some  of  the  basic  elements  of  the  CHA 
technique,  naraaly;1  (1)  a  family  of  target  tracks  are  generated  when  simulated 
and  measured  data  are  cross-correlated;  (2)  one  of  the  tracks  is  the  true 
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target  track  and  the  remainder  are  false  or  alias  tracks.  Although  the  true 
target  track  cannot  be  selected  unambiguously  on  the  basis  of  the  magnitudes 
of  the  correlation  peaks,  the  majority  of  the  alias  tracks  can  be  rejected  by 
ranking  the  correlation  peaks  according  to  size. 


K  *  0  K  *  1  K*  2  K«  3 


Figure  12.  Synthetic  Cross-correlation  coefficient:  tj  =  5.056m,  =  10.0m.  N  =  3.  9 p  -  O.Omr,  r  =  7000.0m. 


Figure  13.  Synthetic  CHA  Target  Tracks:  /;  =  4.8m,  12  -  10.0m,  Op  =  O.Omr.  N  =  7. 
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9.  IDENTIFICATION  OF  THE  TRUE  TARGET  TRACK 

It  has  been  shown  that  family  of  target  tracks  can  be  derived  by 
correlating  measured  data  with  theoretical  data.  The  issue  that  now  needs  to 
be  addressed  is  that  of  identifying  the  true  target  track. 

In  Ref.  1,  Litva  showed  by  means  of  analysis  that  if  the  radar  frequency 
and/or  antenna  height  were  perturbed  in  a  known  manner,  deviations,  exhibiting 
similar  characteristics,  would  appear  on  all  the  tracks  except  for  the  true 
track.  He  did  not  support  his  results  with  experimental  data  but  gave  an 
expression  which  relates  the  magnitude  of  the  deviations  as  a  function  of 
frequency  and/or  height  perturbutions;  namely 


<lw94“{f  +  if +  %  •  -r}' 


(59) 


where  (K-K-j. ) ,  which  takes  on  both  positive  and  negative  values,  identifies 
t'e  order  of  the  track  relative  to  the  true  track  of  number  Rp,  and 


AH  = 


2V  ’ 


Af 

f 


the  relative  change  in  the  radar  frequency. 


Z1 

Az. 


the  relative  change  in  the  antenna  height,  being  its 
average  height. 


3J  1 

”  a  quantity  which  only  becomes  significant  for  very  low 
1  target  heights,  and 


a  curved  earth  correction  factor. 


From  Appendix  A  of  Ref.  1,  an  expression  for  J  in  terms  of  z^,  and  Ag  is 
given  as 


1  -  zA,  f1  -  2»\2]  • 


(60) 


It  follows  from  eqn.  (59)  that  the  sign  of  AHj^  reverses  when  going  from  an 
alias  track  with  (K-Rp)  >  0  to  one  with  (K-Rp)  <  0. 

Equation  (59)  can  be  simplified  by:  (1)  restricting  our  discussion  to 
small  values  of  r,  where  the  flat-earth  approximation  can  be  used  so  that 
J-land  3J/3zp*3zi/j  <<  Azp/zp  and  (2)  keeping  f  constant,  whereby  Af/f«0. 
With  these  approximations, 
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2z. 


r. 


(61) 


If  the  approximations  leading  to  eqn.  (.61)  are  valid,  it  follows  that  for  a 
given  relative  change  in  the  antenna  height  kz\/z\  and  constant  r,  the 
magnitude  of  the  deviations  on  the  Kth  track  is  a  linear  function  of  K. 

From  both  eqn.  (59)  and  its  approximation,  eqn.  (61),  it  follows  that  for 
k-kt,  i AHr| HO. 

Earlier  in  our  discussion  we  saw  that  the  use  of  a  sampled-aperture 
radar  provided  a  built-in  capability  for  perturbing  the  effective  antenna 
height.  This  is  accomplished  by  varying  N  or  dividing  the  array  into 
subarrays.  Therefore,  the  means  for  demonstrating  track  resolution  is  at 
hand. 


9.1  DISCUSSION  OF  SIMULATED  RESULTS 

In  this  section,  simulated  results  are  presented  to  illustrate  the 
antenna-height  perturbation  technique  for  identifying  the  true  track.  The 
sampled-aperture  and  sea-state  parameters  listed  in  Section  8.1  also  apply 
to  the  results  presented  here. 

Figure  14  gives  R[C(r,Z2)]  for  two  antenna  heights  with  r»7000m.  These 
two  heights  were  achieved  by  partitioning  the  array  into  two  sub-arrays  with 
the  uppermost  five  elements  used  for  one  array  and  the  lowermost  five  elements 
for  the  other.  This  gave  effective  antenna  heights  of  z^  *  4.928m  and  z\  - 
4.672m,  where  the  superscripts  U  and  L  denote  uppermost  and  lowermost 
respectively.  With  the  exception  of  the  peak,  K="0,  all  the  peaks  undergo  a 
displacement  when  the  antenna  heights  are  switched  and  the  magnitude  of  this 
displacement  increases  as  K  increases.  A  similar  result  is  obtained  if  only 
the  three  uppermost  (N=3)  and  three  lowermost  elements  are  used.  These 
results  are  given  in  Figure  15  where  in  this  case  z^  and  z^  were  5.056  and 
4.544m  respectively.  "  1 

The  magnitudes  of  the  displacements  in  Figure  15  are  considerably 
larger  than  those  in  Figure  14  because;  (1)  the  relative  effective-height 
variation  is  about  twice  as  large;  and  (2)  the  effect  of  J  in  eqn.  (59). 

The  approximations  leading  to  eqn.  (61)  are  not  accurate  in  describing 
these  results  because  Zj  and  Z2  are  both  small.  Therefore,  the  curved  earth 
correction  parameter  J  must  be  retained.  For  an  accurate  description  of  the 
magnitude  of  the  deviations  in  Figures  14  and  15  one  must  return  to  eqn.  (59). 

If  now  the  results  are  given  in  terms  of  range,  where  target  height  is 
used  in  place  of  the  correlation  coefficient  for  the  dependent  variable,  we 
have  the  results  shown  in  Figures  16  and  17.  These  correspond  to  the  cases 
N=5  and  N=3,  respectively.  Except  for  the  magnitudes  of  the  deviations, 
these  two  results  are  simitar.  Deviations  appear  on  all  tracks  except  for 
the  true  track  (K»0);  their  magnitudes  increase  with  increasing  K  and 
decrease  with  decreasing  r. 


C(r.  Z2)] 


AMK  ■  0  AHK *  1  AHK *  2  AHK ■ 3 


Figure  14.  Synthetic  Cross-correlation  Coefficients;  z^  —  4.928m.  z j  -  4.672m.  N  -  5.  8 p  -  O.Omr, 

r  —  7000.0m. 


Figure  15. 
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As  a  demonstration  of  the  effect  of  J  in  eqn.  (59)  and  the  errors 
incurred  in  setting  J»1  in  eqn.  (61),  let  us  consider  a  particular  example. 
Before  doing  so,  it  should  be  pointed  out  that  one  would  expect  J  to  depart 
from  unity  under  conditions  of  increasing  r  and  decreasing  Z2  simply  because 
the  earth's  curvature  becomes  more  apparent.  In  our  example  z^*4.2m  and 
z->=10m.  When  r=7500m,  J  and  3J/3zj  are  0.805  and  1.376  x  10“^  m-^  respec- 
ively.  For  these  values  the  magnitude  of  3J/3z^*Az^/j  ,'s  only  about  one- 
tenth  that  of  Azj^/z^  in  eqn.  (59)  and  therefore  Justifiably  considered  to  be 
negligible.  The  error  in  | AHjJ  resulting  from  the  approximation  J»1  is  about 
20%.  As  r  decreases,  the  ratio  of  Az^/z^  to  3J/3z^  •  z^/J  increases  and  J 
approaches  more  closely  to  unit.  Therefore,  the  approximations  made  in  eqn. 
(61)  become  more  exact.  For  example,  when  r  is  reduced  to  3500m,  J  becomes 
0.95  and  the  error  incurred  in  setting  J»1  is  reduced  to  5a  and  the  error  in 
neglecting  3J/3zj*Az^/J  becomes  even  less  than  at  r**7500m. 

In  comparing  Figures  16  and  17  one  sees  that  the  ease  of  selecting  the 
true  track  is  enhanced  as  Az^/z^  is  increased.  When  the  CHA  low-angle  filter 
is  used  with  a  radar  for  tracking  low-angle  targets,  it  is  necessary  that 
Az^/z^  or  Af/f  be  chosen  sufficiently  large  so  that  the  track  deviations  are 
detectable  above  the  noise.  For  high  noise  environments,  it  may  be  necessary 
to  implement  matched  filtering  by  correlating  the  noisy  track  -viations  with 
a  replica  of  the  expected  deviations. 


10.  EXPERIMENTAL  ARRANGEMENTS 

The  experimental  data  to  be  presented  in  this  section  are  not  only 
intended  to  demonstrate  the  CHA  technique,  but  also  to  highlight  some 
particular  aspects,  such  as  the  track  Identification  technique  described  in 
the  previous  -ectlon. 

For  this  purpose  a  series  of  experiments  were  conducted  over  an  area  of 
the  Ottawa  River  shown  in  Figure  18  using  an  experimental  sampled-aperture 
array.  The  array  consisted  of  seven  elements,  the  beamwidths  of  which  were 
appro- imately  28°,  with  an  inter-element  spacing  of  d«0.128m  (4A).  The  beam 
maxim;,  of  all  the  elements  were  aligned  in  the  boresight  direction  of  the 
array.  Each  element  was  connected  to  its  own  receiver,  the  local  oscillator 
signal  was  provided  by  a  common  oscillator.  In  this  way  the  phase  relation¬ 
ship  between  the  signals  at  the  antenna  terminals  were  preserved  in  the  I  ano 
Q  signals  appearing  at  the  outputs  of  the  receivers.  Range  information  was 
provided  by  a  Furuno  radar  which  together  with  the  sampled-aperture  array 
formad  the  experimental  radar  facility. 

A  boat  equipped  with  a  mast  which  has  a  CW  source  attached  to  it,  was 
used  to  simulate  low-angle  targets.  The  source,  which  emitted  a  vertically 
polarized  signal  at  9.375  GHz,  was  adjustable  in  height  over  a  range  from  5 
tc  10m  above  the  surface  of  the  water.  An  Alpine  transponder  located  on  the 
same  mast  was  used  in  conjunction  with  the  Furuno  radar  to  provide  accurate 
radar  range  information. 

The  experimental  radar  facility,  which  is  housed  in  a  trailer,  was 
located  at  the  point  shown  in  Figure  18.  The  boresight  of  the  aperture  was 
pointed  down  the  range  towards  Breckenridge.  The  facility  was  under 
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Figure  18.  Site  of  Low-angle  Measurements 


computer  control  whereby  the  radar  range  together  with  phase  and  amplitude 
(1,Q)  information  were  sampled,  digitized  and  recorded  on  magnetic  tape  for 
processing  at  a  later  date.  Further  information  on  the  experimental  arrange¬ 
ments  and  the  characteristics  of  the  sampled-aperture  facility  including  its 
detection  sensitivity  in  the  low-angle  region,  will  be  discussed  in  a  report 
to  be  published  later. 

A  series  of  experiments  were  conducted  whereby  data  were  recorded  for 
various  constant  target  heights  over  a  radar  rar.?e  from  approximately  8000m 
to  5C0m.  From  the  results  of  these  experiments,  ti.ree  have  been  selected  as 
being  representative  of  the  complete  series  and  are  presented  in  the  follow¬ 
ing  sub-section.  However,  before  proceeding  to  this  section,  a  list  of  the 
parameters  which  were  common  throughout  these  experiments  are  as  follows: 

*  4.8m 
SS  -  1 


8 


-12.0  mllx Lradians 
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It  was  intended  that  the  axis  of  the  array  to  b:  coincident  with  the  vertical 
(i.e.,  0p3O),  however,  it  was  found  that  9_  was  decremented  by  12.0  railli- 
radians  relative  to  the  horizontal  line.  This  decrement  was  due  to  a 
combination  of  the  trailer  not  being  level  and  the  sloping  teria<n  over  which 
the  experiments  were  conducted. 


10.1  DISCUSSION  OF  EXPERIMENTAL  RESULTS 

Figures  19  and  20  give  the  amplitudes  recorded  at  two  of  the  array's 
elements  for  a  target  height  of  7.39m,  They  have  not  yet  been  normalized  by 
the  amplitude  recorded  at  a  reference  element.  These  resulcs  illustrate 
classical  interference  patterns  -  in  this  case,  those  appropriate  to  antenna 
heights  5. 18m  and  4.42m  respectively. 

It  is  noted  that  the  period  of  the  interference  pattern  decreases  with 
decreasing  range  in  accordance  with  the  theoretical  results  described 
previously. 

Comparing  Figures  19  and  20  one  observes  a  strong  dependence  of  the 
position  of  the  peaks  and  troughs  of  the  interference  patterns  on  antenna 
height.  This  demonstrates  the  high  sensitivity  obtained  by  using  the 
sampled-aperture  technique  for  deriving  target  heights  by  spatial  sampling 
of  the  interference  patterns. 
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Figure  20.  Measured  Signal  Amplitude  on  Element  7;  Zj  =  4.42m.  if  =  7.39m.  0p  =  -0.012  radians. 

X  =  0.032m. 


The  first  CHA  result  to  be  presented  is  given  in  Figure  21.  It  was 
derived  using  the  same  data  set  from  which  the  examples  in  Figures  19  and  20 
were  taken.  The  measurements,  which  were  taken  at  increments  of  approximately 
2m  in  range,  were  averaged  over  50  data  points  prior  to  processing  in  the  CHA 
filter.  The  data  has  been  subjected  to  some  correction  for  longer  term  drift 
and  instabilities  in  the  sampled-aperture  array.  In  other  words,  it  has  been 
corrected  for  measurement  errors  that  would  normally  be  accounted  for  by 
careful  and  frequent  calibration  of  the  radai  svstem.  The  procedure  used  for 
correcting  the  data  will  be  described  in  more  detail  in  a  subsequent  report. 
The  track  deviations  were  obtained  by  switching,  during  the  course  of  the  CHA 
processing,  from  the  uppermost  three  to  the  lowermost  three  elements.  The 
corresponding  reference  element  heights  were  zr  *  5.06m  and  zr  ■  A. 54m. 


It  should  be  noted  that  the  results  are,  in  general,  similar  to  the 
simulated  results  presented  earlier.  Although  there  are  deviations  on  the 
true  track,  they  are  clearly  smaller  than  those  appearing  on  the  alias 
tracks,  therefore,  the  true  track  is  easily  identified.  The  deviations 
appearing  on  the  true  track  are  probably  due  to  inaccuracies  in  calibration 
and  noise.  Additional  CHA  results  are  given  in  Figures  22  and  23  for  target 
heights  of  3  and  10m.  They  are  similar  to  Figure  21  and  equally  consistent 
with  simulated  results. 

In  the  moael  which  was  used  for  deriving  simulated  data,  the  fact  that 
the  data  were  recorded  over  fresh-water  rather  than  sea-water  was  taken  into 
account  by  using  appropriate  values  of  c/e  and  o  from  Table  1  for  computing 
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the  reflection  coefficient.  Furthermore,  the  array  elerants  were  assumed  to 
be  isotropic  point  receivers.  This  is  a  reasonable  assumption  since,  in  the 
low-angle  region,  the  angle s-of- arrival  of  the  direct  and  indirect  signals 
are  small  compared  with  the  beamwidth  of  the  elements  used  in  the  experimental 
facility. 


Figure  21.  CHA  Tracks  Obtained  From  Measured  Data;  z^  -  5.06m,  z*j  =  4.54m,  Z2  -  7.39m.  N  *  3. 
6p  =  -0.012  Radians,  \  =  0.032m.  <r  =  fl 


Figure  22.  CHA  Tracks  Obtained  From  Measured  Data;  z^  =  5.06m.  z^  =  4.54m,  z^  =  5.0m.  N  =  3, 
0p  =  -0.012  Radians,  A  =  0.032m.  Kf  =  0. 
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Figure  23.  CHA  Tracks  Obtained  From  Measured  Data,  -  5.06m,  Zj  =  4.54m,  -  10.0m.  N 

1 )  =  -0.012  Radians.  X  =  0  032m.  Kj  =  0. 


11.  ALGORITHMIC  PROCEDURE  FOR  TRACK  IDENTIFICATION  AND  DEVELOPMENT 

In  the  previous  section  it  was  demonstrated  that  one  could  pick  out  the 
true  target  track  visually,  simply  by  noting  which  track  had  the  smallest 
deviations.  We  wish  now  to  describe  an  algorithmic  procedure  in  which 
identification  and  subsequent  tracking  of  the  target  may  be  carried  out 
automatically  in  a  digital  computer. 

There  are  two  points  to  note  regarding  this  procedure,  namely;  (1)  we 
must  have  some  knowledge  of  the  magnitude  of  the  deviations  to  be  expected 
at  the  range  of  the  target  and  (2)  we  must  always  maintain  a  good  profile 
of  the  correlation  coefficient  over  the  vector  comprised  of  values  of  z^. 

Implementation  of  track  identification  starts  with  the  first  measure¬ 
ments  of  the  target  range  ^nd  of  the  amplitude  and  phase  distribution  across 
tiie  array.  At  the  onset  estimate  .of  the  deviations  on  the  first  alias 
track,  i.e.,  (K-K^)=l,  is  derived  by  using  eqn.  (61',  which  is 
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The  choice  of  the  initial  group  of  elements  or  subarray  is  unimportant, 
therefore  we  shall  start  with  the  uppermost  group.  Also,  since  wo  have 
little  a  priori  knowledge  at  this  point  of  an  optimum  increment  size  for  z^>, 
we  start  with  a  value  of  lm.  Beginning  at  Z2”lm,  the  complex  correlation 
coefficient  is  calculated  for  each  incremental  value  of  until  22”100m. 

The  values  of  (correlated  heights)  which  correspond  to  peaks  in  the  real 


part  of  the  complex  correlation  coefficient  with  values  that  exceed  0.95  are 
noted  and  recorded.  A  new  increment  size,  *JHC  drived  fr pj  the  average 
separation  in  the  correlated  heights  (*Ay)*  using 
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(63) 


New  height  vectors  for  z^,  corresponding  to  a  new  rmge,  are  derived  in 
preparation  for  the  next  step  in  the  data  processing.  These  values  are 
computed  in  the  vector  range  from  the  IsvfeSt  correlated  height  decremented  by 
zav/4  to  the  highest  incremented  by  z^/4. 

In  the  next  location,  where  the  lowermcst  r,roup  of  elements  are  used, 
both  zAV  and  -me  are  updated  through  multiplying  these  values  by  the  ratio 
of  the  present  range  value  to  that  of  the  initial  range  value.  The  purpose 
of  this  update  is  to  take  into  account  the  decreasing  period  of  the  correla¬ 
tion  coefficient  in  Z2  with  decreasing  range.  In  addition,  jAHK|  is  also 
updated  by  multiplying  it  by  the  same  ratio. 

The  correlation  coefficient  is  then  determined  in  the  manner  described 
previously;  the  height  vectors  for  which  the  correlation  is  a  maximum 
(correlated  heights)  are  noted  and  compared  with  the  previous  ones  in  order 
to  establish  tracks  and  their  respective  accumulated  track  deviations.  The 
previous  correlated  heights  are  then  replaced  with  the  present  ones. 

A  test  is  then  carried  out  so  that  tracks  with  deviations  greater  than 
j  2AHg|  can  be  discarded,  thereby  further  reducing  the  range  of  height  vectors 
in  the  following  iteration.  In  addition,  the  values  of  the  accumulated 
deviations  aie  tested  to  see  if  one  can  be  considered  to  be  a  true  minimum. 

It  is  only  considered  to  be  so  if  its  value  is  smaller,  by  at  least  |AH»|, 
than  anv  of  the  other  accumulated  track  deviations.  The  process  described 
above  is  repeated  until  the  true  track  has  been  identified. 

Once  the  true  target  height,  and  thus  the  true  target  track  has  been 
identified,  a  range  of  height  vectors  ±z.v/4  centred  on  the  true  target 
height  are  then  computed.  The  purpose  of  this  is  to  exclude  the  false  tracks 
during  the  tracking  process. 

In  the  following  iteration,  the  total  number  of  elements  (N-7) 
contained  in  the  aperture  are  then  used  in  determining  the  correlation 
coefficient.  This  has  the  advantage  of  improving  the  signal-to-noise  ratio 
in  the  correlation  process  by  a  factor  which  is  equal  to  the  square  root  of 
the  ratio  of  thp  ..umber  of  elements  currently  being  used  to  that  used 
initially,  thereby  improving  the  accuracy  of  the  estimate  of  the  target 
height.  It  should  be  mentioned  that  once  the  target  track  has  been 
identified,  the  procedure  of  perturbing  the  antenna  height  is  of  little 
further  use. 

In  the  subsequent  iterations  of  the  algorithm,  that  is,  during  the 
process  of  tracking  the  target,  zAy  and  zj^  are  continually  updated.  In 
addition,  a  true  estimate  of  the  target  height  is  determined  by  computing 
a  2-point  running  average  of  the  previous  estimate  ,of  the  target  height  and 
the  present  estimate. 
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In  Che  event  that  the  rack  is  lost  due  to  low  signal-to-noise  ratios, 
the  tracking  process  continues  for  3  iterations  based  on  the  last  known 
estimate  of  the  target  height.  If,  after  this  period,  the  target  track  is  not 
found,  then  the  range  of  height  vectors  is  expanded  to  i*zAV,  centered  on  the 
last  known  estimate,  and  the  procedure  of  perturbing  the  antenna  height  is 
reintroduced  so  that  the  true  target  track  can  be  re-identif ied . 


11.1  RESULTS 

The  final  stage  of  CUA  processing  consists  of  track  selection.  A 
demonstration  of  track-selection  is  given  in  Figures  24  to  26  using  the  sane 
data  that  was  used  in  Figures  21  to  23.  It  should  be  noted  that  in  each 
instance  the  algorithm  has  identified  the  correct  target  track  and  has 
successfully  developed  the  track  over  the  entire  range  for  which  data  was 
available.  Bearing  in  mind  that  the  nominal  target  heights  were  respectively 
7.39,  3.0  and  10.0  metres,  it  is  readily  seen  that  the  tracking  accuracy  of 
the  algorithm  is  very  high  indeed. 


12.  SUMMARY  AND  CONCLUSIONS 

In  this  report  is  given  a  review  of  the  CHA  low-angle  tracking 
technique.  In  so  doing,  more  analytical  results  are  presented,  thereby  more 
clearly  demonstrating  the  basis  for  the  algorithm  and  documenting  the  various 
stages  of  processing.  One  stage  of  the  processing  that  receives  considerable 
attention,  because  of  recent  developments,  is  the  track-selection  algorithm. 
This  algorithm  is  based  on  a  perturbation  technique,  where  either  the  antenna 
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Figure  24.  CHA  Target  Track  Obtained  From  Measured  Data:  =  7.39m 
(Using  the  Same  Data  as  m  Figure  21). 
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Figure  25.  CHA  Target  Track  Obtained  From  Measured  Data:  i2  =  5  0m 
(Using  the  Same  Dicta  as  in  Figure  22). 


Figure  26.  CHA  Target  Track  Obtained  From  Measured  Data;  i2  =  0.0m 
(Using  the  Same  Data  as  in  Figure  23). 


height  or  RF  frequency  is  purposely  varied.  It  is  shown  that  a  change  in 
antenna  height  of  die  order  of  11%  was  sufficient  to  permit  successful 
operation  of  this  algorithm.  Basically,  this  antenna-height  variation  is 
sufficiently  large  to  produce  deviations  on  the  false  tracks  which  are 
larger  than  the  deviations  on  the  true  track  caused  by  errors  in  calibration 
and  noise. 

A  considerable  amount  of  the  discussion  is  centred  on  the  characteris¬ 
tics  of  the  interference  pattern  at  the  radar  resulting  from  the  direct  and 
indirect  rays  from  the  target.  This  subject  was  chosen  as  a  topic  of 
discussion  so  as  to  set  the  stage  for  an  introduction  to  sampled-aperture 
technology,  in  particular,  the  use  of  a  sampled  aperture  array  for  recording 
data  which  can  subsequently  be  used  for  testing  and  developing  the  CHA 
algorithm.  The  algorithm  as  described  here  is  recast  somewhat  from  that 
described  by  Litva1  so  as  to  be  compatible  with  data  recorded  with  an 
experimental  sampled-aperture  radar  facility.  This  report  contains  a 
thorough  description  of  the  analytical  processes  inherent  to  the  revised 
algorithm  and  its  software  implementation. 

Finally,  the  algorithm  is  tested  with  experimental  data  recorded  over 
the  Ottawa  River  during  che  summer  of  1980.  It  is  shown  that  with  good 
calibration  of  the  sampled-aperture  facility,  the  CHA  algorithm  is  capable 
of  tracking  low-altitude  car gets  for  ranges  between  7000  and  500m.  In 
addition,  rms  tracking  accuracies  of  2.0m  or  0.5  milliradians  are  demonstrated 
with  target  altitudes  at  least  as  low  as  5.0m. 
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